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Abstract. We investigate the localization properties of a one-dimensional Dichro- 
matic optical lattice in the tight binding regime, by discussing how exponentially lo- 
calized states emerge upon changing the degree of commensurability. We also review 
the mapping onto the discrete Aubry-Andre model, and provide evidences on how the 
momentum distribution gets modified in the crossover from extended to exponentially 
localized states. This analysis is relevant to the recent experiment on Anderson local- 
ization of a noninteracting Bose-Einstcin condensate in a quasiperiodic optical lattice 
[G. Roati et ai, Nature 453, 895 (2008)]. 
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1. Introduction 

One-dimensional quasiperiodic potentials have been the object of an extensive 
theoretical investigation since the 80's [UE1I3JIHI51E1II1E1I91 E2; owing to their 
interesting feature of having a spatial ordering that is intermediate between periodicity 
and disorder. The paradigm of this class of models is the so called Harper pQ or Aubry- 
Andre model [2], described by the following hamiltonian 



that, for a certain class of irrational values of/? [51 E], displays a "metal-insulator" phase 
transition from extended to exponentially localized states, just as Anderson localization 
in random systems in higher dimensions [21 HH [12]. This model has the interesting 
properties of being self-dual under a suitable transformation to momentum space, the 
self-dual point being A/ J = 2 at which the transition occurs [2] . 

This kind of hamiltonians can now be engineered with ultracold atoms loaded 
in optical lattices [131 [EH H51 EE CEO EES], opening the way to study the localization 
properties of quasiperiodic systems from the experimental point of view. The first 
experimental observation of Anderson localization in quasiperiodic potentials has been 
recently reported in [T9l 120] . This experiment has been realized by using a non 
interacting Bose-Einstein condensate loaded in a bichromatic potential, obtained by 
superimposing two one-dimensional optical lattices with different wavelengths. 

In this paper we investigate the localization properties of the bichromatic potential 
employed in current experiments, discussing the relationship with the discrete model in 
([TJ, and also the effect of the finite size of the system due to the presence of an additional 
harmonic confinement. We will provide evidences on how exponentially localized states 
emerge upon increasing the degree of incommensurability of the system, and how the 
transition from extended to localized states can be extracted from the momentum 
distribution, that is easily accessible in current experiments from the absorption images 
after a time of flight [19] . 

2. The bichromatic lattice 

Recent experiments with cold atomic samples feature a one-dimensional bichromatic 
potential obtained by superimposing two optical lattices, of the form [161 EE] 



where fcj = 2ir/Xi (i = 1,2) are the lattice wavenumbers, Sj are the heights of the two 
lattices in units of their recoil energies Em = h 2 / (2mAf), and is an arbitrary phase. In 
general, the potential of wavelength Ai is used to create a tight binding primary lattice, 
of period d = Ai/2, that is weakly perturbed by the secondary lattice of wavelength A 2 . 




(1) 




(2) 
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The single particle properties of the potential Vb(x) are obtained by diagonalizing 
the hamiltonian 

H=-^Vl + V b (x), (3) 

that can be conveniently rewritten in units of main lattice parameters (fcf 1 and Em as 
units of lenghts and energies, respectively) as 

Hj E m = - V 2 + Sl sin 2 (0 + s 2 p 2 sin 2 (/3£ + 0) 

= H 1 + s 2 p 2 sin 2 (P^ + <P) (4) 

with £ = kix, (3 = Ai/A 2 and E R2 /Em = P 2 - 

It can be easily demonstrated that the presence of the second lattice does not change 
substantially the position of the minima but instead shifts the energies in a range of size 

s 2 (3 2 02! • 

2.1. The tight binding limit: mapping onto the Aubry- Andre model 

In the tight-binding limit this system can be mapped to the Aubry-Andre model [2] 
by expanding the particle wavefunction ip over a set of Wannier states, maximally 
localized at the minima of the primary lattice (from now on "lattice sites", labelled 
by j) Wj = w(£ — £j), in the lowest band, = J2j c j\ w j)- The Aubry-Andre model can 
be written in terms of the projector over the Wannier states as 

H — > \wi)(wi\H\wj){wj\ (5) 

ij 

or in terms of the Cj coefficients 

H — = ^c*c>,|#K-)- (6) 

ij 

Let us evaluate explicitly the expectation value of the hamiltonian on the Wannier 
basis; in the tight binding limit we can make the following approximation 

(wi\H\ Wj ) w eoSy - J5 i>j±1 + 5 ir s 2 [3 2 J d£ sm 2 (p£ + 0)K(OT (7) 

with 

e = J d^ Wi (0H lWi (0 ; J = J di w i+1 (OH lWi (0 (8) 

where we have retained only the onsite contribution of the secondary lattice [T71 [18] 
and the tunneling between nearest neighbouring sites [21 J (for shallow lattices, the 
interaction between Wannier states beyond nearest neighbours becomes relevant too 
[10J). The integral in last term of Eq. ^ can be rewritten by using the relation 
sin 2 (/?£ + <p) = (1 - cos(2/?£ + 2<f>))/2, and noticing that 



d£ cos(2/3£ + 0')K(£)| 2 = cos(2vr/5i + 0') / d£ 



cos(2/30k(0| 2 - (9) 
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Finally, by dropping constant terms and defining 

A S ^ 



dC cos(2/30KOr ( 10 ) 



we have 



(wi\H\wj) ~ —J5ij±i — 5ijAcos(27ij3i + <fi') (11) 

Therefore, in the tight-binding approximation the hamiltonian for the bichromatic 
potential can be cast in the form of the Aubry- Andre model [21 13] 

H = -jJ2 [ c )+i c 3 + c* c i+0 + A J2 cos ( 27r A? + ^) c i c i ( 12 ) 
i j 

or 

if = - jJ2(\wj){w j+1 \ + \w j+1 )( Wj \) + A^cos(27r#7 + 0')K')MM) 
where the characteristic parameters J and A are connected to the parameters directly 



tunable in the experiment, s±, S2 and f3, via Eqs. (|8j), (10). These can be 
estimate by using a Gaussian approximation for the Wannier functions, |«;(£)| 2 ~ 
(l/ v /vr)si /4 exp(- v /sl^ 2 ), that gives J = (A/^/tt)s^ 75 exp(-2 v /sl) [22J and A = 
(s2/3 2 /2) exp(— p 2 /y/s~i) [ID] . A more precise expression can be obtained by calculating 



numerically the integrals in Eqs. (|8j), (10) in terms of the maximally localized Wannier 
functions. This yields [23] 

J ~ 1.43s?- 98 exp(-2.07v^) , (14) 



while for A we can make the following ansatz 

A * 'jfe-PW (15) 

with a and 7 to be determined from a fit of the numerical evaluation of the integral 
in (10). We note that in the literature it is common to find A defined without the 
exponential factor [181 EE US]; though this is a legitimate definition of the intensity of 
the secondary lattice, the contribution of the exponential correction may be necessary 
for a precise mapping between the continuous and the discrete models, depending on 
the values of s\ and (3. 



2.2. Tuning the degree of commensurability 

A common choice in the study of the Aubry-Andre model is the inverse of the golden 
ratio, (3 = tp^ 1 = (Vo — l)/2, for which the model displays a "metal-insulator" 
phase transition from extended to localized states at A/J = 2 [2]. In particular, the 
insulator phase is characterized by the absence of mobility edges, all the eigenstates 
being exponentially localized with the same localization length I = l/ln(A/2J). This 
case has been extensively investigated in the literature, see e.g. [31 El El [S] • 

From the mathematical point of view, a necessary and sufficient condition to display 
a sharp transition is j3 to be an irrational Diophantine number [B]. However, in the 
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Figure 1. Plot of the bichromatic potential for s x = 10 (J « 0.02), A/J= 10. The 
blue stripes indicate the first and second Bloch bands of the primary lattice, whereas 
the red one represents the amplitude 2A of the perturbation induced by the secondary 
lattice. Notice that the latter is on an energy scale much lower that the bottom of the 
lowest band. 



real world one has to be aware of two facts: i) real system have a finite size; ii) 
experimentally (and numerically) the wavelenghts can be given with a finite number 
of digits (in practice all number are rational). This means that what really matters 
is not the distinction between commensurable and incommensurable, but the actual 
degree of commensurability and the ratio between the periodicity of the potential and 
the system size (fixed by the boundary conditions). In practice, in order to observe a 
transition from extended to localized states it is sufficient to have a large enough number 
of lattice sites within the actual periodicity of the potential, and a system size that does 
not greatly exceed the latter (in order to avoid periodic replicas). 

To illustrate how the properties of the system change by varying the degree of 
commensurability, we consider the solution of the full hamiltonian in Eq. (|4]), obtained 
by solving the Schrodinger equation (for convenience here we fix = 0) 

+ Sl sin 2 (0 + s 2 (3 2 sm 2 ((30] V>(0 = E ^(0- (16) 

Here we focus on the tight-binding regime, and throughout this paper we fix s\ = 10 
(J ~ 0.02). With this choice, an accurate expression for A is given by 

A * 'J^e-aF (17) 

where the parameters a and y can be obtained by a numerical fit, yielding a = 0.385, 
y = 1.95. In Fig. [l]we plot a sketch of the potential for A/J = 10 (corresponding to 
s 2 ~ 0.47), compared with the first and second Bloch band of the primary lattice. 
This picture shows that the large separation between the first and second bands 
indeed justifies the mapping onto the single band Aubry-Andre model, and that the 
perturbation introduced by the secondary lattice is at an energy scale much below the 
lowest band. 
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Figure 2. Color plot of the inverse participation ratio (-P -1 ) for the Dichromatic 
potential as a function of A/ J and the ratio (3 of the two wavelenghts. Black regions 
correspond to extended states (P^ 1 ~ 0), whereas those colored in orange/yellow to 
localized states (P _1 > 0.5). The size of the intermediate palette of colours measures 
the width of the crossover. The labels indicate the position of simple ratios of the two 
wavelenghts. 



Differently from the common choice of periodic boundary conditions [3j |9], here 
we consider instead vanishing boundary conditions to account for finite size effects, in 
view of the fact that in the experiments there is always a residual harmonic confinement 
that introduces a natural cutoff at large distances (this choice does not affect the bulk 
properties). The specific effect of the harmonic confinement will be discussed later on. 

To illustrate how the properties of the system change by varying /3, in Fig. [2] 
we plot the inverse participation ratio = J d£\i[)(£)\ , that measures the inverse 
of the number of occupied lattice sites [8]. This picture shows that it is possible to 
have localization for a large enough A/ J, in a wide range of values of j3, and that the 
transition broadens in the neighbourhood of simple rational values. 

2.3. A Fibonacci sequence 

To get further insight on how the localization transition changes when we increase the 
incommensurability of the system here we consider a set of bichromatic potentials whose 
wavelengths are in the ratio f3 n = f n / fn+i of two consecutive numbers of the Fibonacci 
sequence, defined by the recursion relation f n+ i = f n + fn-i with fo = fi = 1 (each 
number after the first two is the sum of the previous two numbers), /„ =1, 1, 2, 3, 5, 
8, 13, 21, 34, 55, 89, 144, 233, etc. [5]. As an example here we choose A/J = 10, for 
which we expect localization in the incommensurate limit (3 = ip~ l [8]. 

In Fig. [3^, we plot the groundstate density of the bichromatic potential, for 
Po = 1/1, p 4 = 5/8, (3 8 = 34/55, (3 12 = 233/377, for a system of length L = 500 (number 
of lattice sites) . This picture shows how the density distribution of the groundstate 
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Figure 3. Plot the groundstate density for the Dichromatic (a) and random (b, see 
j pl| potential for O = 1/1, /3 4 = 5/8, /3 8 = 34/55, f3 12 = 233/377 (from top to 
bottom), and A/ J — 10, for a system of length L = 500 (number of lattice sites). 



gets modified as we increase the actual periodicity of the system, with peaks localized 
every f n+ \ sites of the main lattice, characterized by an exponential decay away from the 
localization center over an increasing number of lattice sites, ^(Ol 2 ~ ex P( — 1£ — £o|/0- 
As the ratio f3 of the two wavelengths approaches an irrational number the distance 
between two consecutive localized states goes to infinity, and in the limit only one 
state survives. Note that the localization is not due to the suppression of tunneling 
- for (3 rational the states are extended - but to the randomization of onsite energies 
characteristic of Anderson localization [3]. 



2.4- Quasiperiodic vs random 

To compare the behaviour of the quasiperiodic lattice with that of a pure random system 
we consider a set of disordered lattice potentials, but with a periodicity of f n+ \ sites 
as before, obtained by adding to the primary lattice a sequence of f n+ i on-site energies 
€j extracted from a uniform random distribution in the range[— A, A]. This system is 
described by the tight binding hamiltonian 

H = - J Yl (K)H+il + K+i)K'l) + ^2 e j\ w j)( w j\ ( 18 ) 



J J 



obtained by replacing the second term in Eq. (12) by a random on-site energy ternjf] 
In the limit n —>■ oo this system realizes the one-dimensional Anderson model [TT] with 
nearest-neighbour coupling. The corresponding groundstate distributions for A/J = 10 

| In practice we have solved the continuum version of this model, obtained by shifting the minima of 
the primary lattice by tj. 



Exponential localization in one-dimensional quasiperiodic optical lattices 



8 



are shown in Fig. [HJo. This picture evidences a very similar behaviour of the two 
models, suggesting that the mechanism of localization in quasiperiodic potentials is 
essentially the same of Anderson localization in disordered systems, irrespective of 
the fact that the shuffling of the energies among different lattice sites occurs in a 
correlated or uncorrelated way (for the quasi-periodic and random cases, respectively) 
[I]. This correspondence is remarkable, and supports the use of bichromatic potentials 
as effective tool for implementing quasi-disorder in current experiments with ultracold 
atoms [131 EE HE! EE]. The main difference is that for pure random system, in one 
(and two) dimension all states are always localized, for any amount of disorder [8j [12] . 
Instead, one-dimensional quasiperiodic potentials may show both extended and localized 
states, displaying a "metal-insulator" transition as a function of the strength of disorder 

0H- 



3. Finite size systems: the effect of harmonic confinement 

The localization properties of a bichromatic potential have been recently demonstrated 
experimentally in [19], where it has been reported the first observation of Anderson 
localization for a noninteracting matter wave. In this section we analyze the localization 
properties of this specific system, discussing the broadening of the "metal-insulator" 
transition for the specific choice of the wavelenghts, the effect of the harmonic 
confinement, and the modification of the momentum distribution across the transition. 

In the experiment (19] a noninteracting Bose-Einstein condensate of 39 K atoms is 
initially prepared in a 3D harmonic potential Vh {x, r±) = muj 2 x 2 /2 + mu> 2 _r 2 _/2, r± 
being the radial coordinate in the transverse direction (with respect to the lattice), and 
then loaded in a bichromatic potential of wavelengths Ai = 1032 and A2 = 862, yielding 
(3 = 1.1972 .... In the most commensurate hipotesys (neglecting decimal digits in the 
wavelenghts), this choice corresponds to a strict periodicity every n\ = 432 sites of the 
primary lattice. Since the condensate is non-interacting, the problem is separable and 
the axial and radial wavefunctions factorize. In this case the presence of the bichromatic 
lattice affects only the axial wavefunction, that can be obtained by diagonalizing the 
hamiltonian 

ft 2 1 

H = -— V 2 x + V b {x) + -mu 2 x 2 . (19) 
2m 2 

The only difference with what discussed in the previous sections is the presence of 
the axial harmonic confinement!!] , that fixes the size of the system to a value of the order 
of the oscillator length, a^o = yh/muj. In the following we will assume dho *C in 
order to consider the case of systems whose size does not exceed the overall periodicity 
of the bichromatic potential. Moreover, we may expect the localization properties to be 
not substantially affected by the harmonic trap when the system extends over several 

§ In the numerical calculations we have also included a tiny phase shift between the bichromatic and 
harmonic potentials in order to avoid perfectly symmetric solutions, that may be quite rare in the 
experiment. 
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Figure 4. Left: participation ratio of the groundstate of the bichromatic potential of 
Ref. [19] as a function of A/J, for a trapping of frequency u res = 2ir x 0.5 Hz (red, 
continuous line) and u> = 2ir x 5 Hz (red, dashed line). The dotted line represents the 
participation ratio for (3 — (\/5 — l)/2, in the case u> = ui res . Right: density plot of 
the groundstate for oj = u) res . 



lattice sites, that is when aho 3> d. In this case the harmonic trap only affects the 
boundary conditions but not the localization behaviour in the bulk. In the experiment 
the lowest value of the trapping frequency is fixed by the residual trapping due to the 
focusing of the laser beams that creates the primary lattice, uj res ~ 2tt x 0.5 Hz. Instead, 
the measurements of the momentum distribution have been taken with uj rs 2tt x 5 Hz. 
Both these values satisfy the above requirements, and correspond to a^ Q / d = 44, 14 
respectively. 

In Fig. [4] we show the participation ratio and the density plot of the groundstate 
of the bichromatic potential as a function of A/ J, for the residual trapping, uj = uj res . 
We note that in [19] the intensity of disorder has been defined according to the common 
choice A = S2(3 2 /2 [181 HE]; here we will instead use the definition in Eq. (17) for a 
more precise comparison with the discrete Aubry- Andre model. This picture shows that 
for the particular choice of the wavelenghts of [TH] there is a broadening of the "metal- 
insulator" transition with respect to the "incommensurate" case (5 = <fi~ l , and that the 
eigenstate size becomes comparable with the lattice spacing at about A/J = 5. 

The modification of the spectrum as a function of A/J is shown in Fig. [5^b 
for both values of the trapping frequency. Here we have considered the energy range 
corresponding to the lowest band of the primary lattice. This picture illustrates that, in 
case of weak confinement (uj = 2tt x 0.5 Hz), several mini gaps open in the spectrum as 
we turn on the secondary lattice, the main ones being those corresponding to the beating 
of the two lattices [15] . Instead, for a tighter trapping (uj = 2ir x 5 Hz) the spectrum is 
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Figure 5. Top row: spectrum of the Dichromatic potential as a function of A/J, for 
u> = 2ir x 0.5 Hz (a) and uo — 2tt x 5 Hz (b). (c),(d): the corresponding ground states 
(at the center, red continuous line) and the first six excited states (dashed lines) for 
A/J = 10. 



deeply affected by the presence of the harmonic potential, that produces an enlargement 
of the band size and a different redistribution of the energies. In particular, the presence 
of the harmonic confinement affects the spatial arrangement of the localized eigenstates, 
as shown in Fig. [5jxl. For a shallow enough confinement the localized eigenstates are 
shuffled all along the system, whereas a tighter confinement produces a hierarchical 
arrangement starting from the trap center, with eigenstates neighboring in energy that 
are disposed every five sites (with a reflection symmetry about the trap center). This 
feature allows to tune the arrangement of the eigenstates around the trap center by 
changing the confinement frequency, without affecting their exponential localization 
behavior [19] . 

3.1. Momentum distribution 

Important indications about the modification induced by the secondary lattice on 
the eigenstates of the system can be obtained from their momentum distribution. 
This is shown in Fig. [6] where we plot the groundstate density |VK£)| 2 (left) and 
its momentum distribution p(k) = \ifj(k)\ 2 (right) for increasing values of A/J, with 
V'(^) = / d^,ip(C)e ■ Here we consider the same value of the trapping frequency 
uj = 2tt x 5 Hz used in [19J to make the analysis of the momentum distribution after the 
time-of-flight. 
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When only the primary lattice is present, p(k) displays the typical interference 
pattern with three peaks at k — 0, ±2&i reflecting the periodicity of the system and the 
fact the the eigenstates are extended. The tiny width of the peaks indicates that the 
wavefunction spreads over many lattice sites [21] . When the secondary lattice is added, 
additional momentum peaks appear, corresponding to the beating of the two lattices. 
The presence of these peaks can be easily explained for a commensurate system defined 
by ni\i = n 2 X 2 (here we consider n 2 > ni = L as in the experiment, (3 = r^/rii)) by 
transforming the discrete model in (13) in the momentum representation (ki = 21 /L in 
units of ki) 

L 

H = -2JJ2 C0S { 2 * 1 / L )\ k i)(ki\ + Aj2$i-i>\ki)(ki>\ (20) 

1=1 w 
where we have defined \h) — (1/y/L) J2j G mkl ^\wj), and 

L L 

$ lv = _ ^ ^2 e~ t27Tj{{l - l ' )/L±l3) = —^2^2 e~ i7r{k ' kl ' ±2k2)j (21) 

2L 3=1 ± 2L j=l ± 

with k 2 = (3ki = j3. It is easy to prove that the term <5;_//, that determines the coupling 
between different momentum components, is non vanishing when \ki — ki>\ = 2k 2 — 2k\ 
or \ki — ki'\ = 4&i — 2k 2 . When going back to the continuous model, the above relations 
correspond to the additional peaks in the momentum distribution, that appear at a 
distance ±2(k 2 — k\) around those at integer multiples of ±2k\ due to the primary 
lattice. 

By further increasing the intensity of the secondary lattice, the momentum 
distribution p{k) broadens and its width eventually becomes comparable with that of the 
first Brillouin zone. This take place in correspondence of the appearance of exponentially 
localized states whose extension shrinks below that of a single lattice site. The residual 
modulation of the momentum distribution in the two bottom profiles of Fig. [6] (right) 
indicates that the localization is nontrivial, in the sense that the tails of the eigenstates 
extend over several lattice sites even for large A/ J (as shown in the left column). This 
analysis accounts for the behavior observed in [Tj5] across the localization transition. 



4. Conclusions 



We have analyzed the single particle localization properties of a one-dimensional 
bichromatic potential obtained by superimposing two optical lattices of different 
wavelenghts, discussing how the degree of commensurability of the two lattices affects 
the appearance of exponentially localized states, in the tight binding regime. We have 
shown that in order to observe a transition from extended to localized states in a finite 
size system what really matters is not the distinction between commensurable and 
incommensurable, but the possibility to have a large enough number of lattice sites 
within the actual periodicity of the potential. 
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Figure 6. Plot of the groundstate density (left, in log scale) and of its momentum 
distribution (right, note the different scale of the vertical axis) for different amplitudes 
of the secondary lattice, A/ J — 0, 2, 5, 10, 20 from top to bottom, in the presence of a 
harmonic trapping of frequency u> — 2tt x 5 Hz. 



We have also reviewed the mapping onto the discrete Aubry-Andre model [2] and 
made a comparison with pure random disorder, providing evidence of the similarity of 
the localization mechanism in quasiperiodic and random systems, in support of the use 
of bichromatic potentials as effective tool for implementing quasi-disorder in current 
experiments with ultracold atoms [131 E51 ESI EH] ■ 

Finally, we have discussed the effect of an additional harmonic potential and 
analyzed the modification of the momentum distribution across the "metal-insulator" 
transition, making a direct discussion of the experiment [TH], where it has been 
reported the first observation of Anderson localization of a noninteracting Bose-Einstein 
condensate in a quasiperiodic optical lattice. 
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